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ABSTRACT 

The orbits of stars close to a massive black hole are nearly Keplerian ellipses. Such orbits exert 
long term torques on each other, which lead to an enhanced angular momentum relaxation 
known as resonant relaxation. Under certain conditions, this process can modify the angular 
momentum distribution and affect the interaction rates of the stars with the massive black 
hole more efficiently than non-resonant relaxation. The torque on an orbit exerted by the 
cluster depends on the eccentricity of the orbit. In this paper, we calculate this dependence 
and determine the resonant relaxation timescale as a function of eccentricity. In particular, we 
show that the component of the torque that changes the magnitude of the angular momentum is 
linearly proportional to eccentricity, so resonant relaxation is much more efficient on eccentric 
orbits than on circular orbits. 
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1 INTRODUCTION 

It is now commonly accepted that many galaxie s contain mas- 
sive black holes ( MBHs) at their centres (e.g., IGebhardt et al.l 
l2003l : lMilleill2006h . with masses 10*' Mq > M. > IQ^Mq- The 
proximity of our own galactic centre, at a distance of 7.62 ± 
0.32 kpc teisenhauer et alfcoOSi) . allows for the astrometric study 
of i ndividual stellar orbits in the potential well of the MBH 
(e.g.. lSch6del et al.l2002l ; lGhez et al.l2"003l ; teisenhauer et alj2005h . 
These observations showed that the galactic MBH mass is M, = 
(3.61 ± 0.32) X 10*^ Mq teisenhauer etai]|2005h . The presence 
of very dense stellar clusters near MBHs leads to a wide variety 
of interesting phenomena. For a compr ehensive re v iew o f stellar 
phenomena in the Galactic centre, see I Alexander! ( |2005|) . Many 
of these phenomena arise because the two-body relaxation time is 
shorter than the age of the systems, so that the orbits of stars are 
significantly redistributed during their lifetimes. For example, the 
relaxation time in the Galactic centre is estima ted to be several Gyr 
( lAlexandeill 19991 : [Alexander & Hopmaij|2003h . 

A common assumption in stellar dynamics is that the mech- 
anism through which stars exchange angular momentum and en- 
ergy is dominated b y uncorrelated two-body interactions (e.g., 
IChandrasekhaj 119431) . The orbits of the stars are largely deter- 
mined by the potential of the smoothed density and the devi- 
ations from this potential caused by the individual stars lead 
to perturbations that evolve the orbits (Henon 1973). Systems 
with MBHs at their centres cannot (yet) be studied by di- 
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rect jV-body integrations using a realistic number of particles. 
Hence the methods used for studying the dynamics of these 
systems rely on the above assumptions. Such st udies include 
Mont e Carlo simulations (e.g .,[shapiro & Marchan^ |T978: Freita3 
1200 ll : iFreitag & Bend I2OO2I : iFreitag 2003; Freitag et al. 200^ 
and Fokker-Planck methods (e.g., [Sahcall & Wolf ,1976i. 1 19771; 
lMurDhvetai]|l99ll : iHopman & Alexanderll2006bh . 

Near a MBH, the potential is nearly Keplerian, which leads 
to closed elliptic orbits. There are two main reasons for deviations 
from Keplerian orbits: the contribution to the potential from the 
stars and general relativistic effects. However, for most orbits where 
the potential is dominated by the MBH, the time scale for preces- 
sion is large, and the orbits remain nearly stationary with respect to 
each other over many periods. Consequently, the assumption that 
the intera ction between the stars ar e independent is not valid, as 
shown bv' Rauch & Tremaind jl99^. They argued that a better de- 
scription is given when the interactions are considered to be be- 
tween different orbits, rather than between different point particles. 
Averaged over many periods, one can think of the mass of a star 
being smoothly distributed over the orbit, with the linear density 
in a small segment proportional to the time the star spends in this 
segme nt. The orbits then form massive "wires". Ranch & Tremain^ 
( Il996h showed that the torques between the wires lead, under some 
circumstances, to very efficient angular momentum relaxation. 

Part of the interest in efficient angular momentum relaxation, 
is due to interesting phenomena which occur when stars have very 
close interactions with MBHs when they are on highly eccentric 
orbits. Such phenomena include the tidal disrup t ion of stars (e.g., 
iFrank & Reeslll976l : iLightman & Shapi"rolll977l ; lReeslll988l) . and 
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the e mission of gravitational waves with detectable frequencies 
(e.g..lSigurdssonlll99l lFreitag|l2003l : iHopman & AlexandeJiiool 
l2006air Since resonant relaxation increases the rate of angular mo- 
mentum scattering, star s reach highly eccentric orbits more rapidly. 
iRauch & Ingailsl ( Il998l) studied the consequences for the tidal dis- 
ruption rate, while iHopman & Alexander! i2006ah considered the 
enhancement of the rate at which gravitational wave sources spiral 
in. Resonant relaxation appears to be more important for the latter 
case, since it occurs closer to the MBH, where deviations from a 
1/r potential are smaller. 

Resonant relaxation affects both the magnitude and the direc- 
tion of angular momenta. For processes that require orbits of high 
eccentricity, only the part which af fects the overall ma gnitude of 
the angular momentum is of interest. lRauch & Tremainel ( 11996) es- 
timate the torques of a stellar cluster on a test star, and deduce the 
relaxation rate from this. They do not consider the dependence of 
this process on the eccentricity. This is potentially important, as can 
be seen from the following example. 

Consider a star in a circular orbit in the sy-plane, and a mass 
at a general point p= (x,y,z).\n this case, the angular momentum 
of the star is in the z direction, J = Jz. In the approximation that 
the orbit does not change, the torque r that the mass at p exerts on 
the star would have no z component since the contributions from 
one half of the orbit will cancel the contributions from the other 
half, T = {Tx,Ty,Q). Since r-J = 0, the mass at p can rotate 
the orbit, but never affect the magnitude of the angular momentum, 
and since p was a general point, this is true for all other points as 
well. We can therefore conclude that resonant relaxation can never 
modify the eccentricity of an e = orbit. 

For any orbit with e > 0, resonant relaxation will change the 
eccentricity. In this paper, we calculate how the efficiency of reso- 
nant relaxation depends on the eccentricity of the orbit of a given 
star. 



For a large number of stars in the region nea^j a, the sum 
of the torques will nearly average out to zero. However, because of 
statistical fluctuations, there will be an excess torque in an unknown 
direction of order rjv ~ VN ti . If the orbit of the test star lies in 
the xy-plane, so that its angular momentum is in the z-direction, 
only the z-component of the torque can affect the magnitude of the 
angular momentum (or the eccentricity) of the star. The resulting 
an gular momentum ch anges were called scalar resonant relaxation 
bv lRauch & Tremainei (1996). If the eccentricity of a test stars is e, 
then the typical net torque in the z-direction will be 



r,(a,e)=Ps{e)VN 



(3) 



where f 3s(e) is a dimension l ess fu nction of eccentricity. 

In iRauch & Tremain^ ( 1 19961) and later papers, the ecce ntric- 
ity dependence of /3s was ignored. iRauch & Tremain^ ( Il996h per- 
formed A''-body simulations to find the eccentricity averaged value 
of l3s{e), and found that for an isotropic cluster with eccentricity 
distribution 



Mso(e) = 2ede, 



(4) 



the ty pical value of (3s is (3a — 0.53 ± 0.06 jRauch & Tremairi^ 
1 19961 table 4d). However, from the discussion in the introduction, 
it follows that /3s (0) =0. 

After a time t^,, the orientation of the test star with respect 
to the ambient cluster changes. This can happen either because its 
orbit has precessed, or because the orbits of most of the stars around 
it have precessed. Over this time, the angular momentum would 
change by 



/3s(e)%/iV^t^(a,e) 



(5) 



2 RESONANT RELAXATION 

In a Keplerian orbit there is a 1 : 1 resonance between the angu- 
lar and radial frequency and hence the orientation of the orbit is 
fixed in space. Because of the general relativistic effects, all orbits 
around a MBH actually precess and are never exactly Keplerian; 
in addition, the contribution to the potential from the star cluster 
around the MBH will also lead to precession. However, for nearly 
Keplerian systems, the orientation of the orbit with respect to other 
orbits can be assumed to be fixed in space over some time t^^ 3> P, 
where 



P{a) = 2tt 



\GM. 



1/2 



(1) 



is the period of the star. 

Over a ti mescale P <^t ^t^j, the st ars can be represented as 
massive wires l lRauch & Tremainei [l996l) with the mass smeared 
out over their orbits. These wires exert mutual torques on each 
other. The magnitude of the torque exerted by a star of mass 
and semi-major axis a on another star with equal semi-major axis 
is estimated byQ 



(2) 



For t > tui the torques on a particular star-wire become ran- 
dom, and the change in angular momentum grows as a random 
walk. The resonant relaxation time Trr is defined as the time 
is takes for a star to have its angular momentum changed by an 
amount of the circular angular momentum 



Jc(a) = VGM.a. 



(6) 



Since it takes (Jc/AJi^)^ random steps to make this excursion in 
angular momentum space, and each step takes a time t^j, the reso- 
nant relaxation time is given by 



Trr (a, e) 



Jc 



27r/3s(e) 



M.y 1 P(a)^ 



(7) 



In the following section we describe a method to find the 
torque on a wire which we use to determine the eccentricity de- 
pendence of /3s . We note that from equation Q is can be seen that 
since /3s (e) — ^ for e — > 0, scalar resonant relaxation becomes 
very ineffective for nearly circular orbits. 



^ For brevity, we use angular momenta, torque etc. per unit mass. 



^ In Section[3]we show that the maximal distance at which stars still have 
a large contribution to the torque is 2a. 
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3 THE WIRE APPROXIMATION FOR TORQUE 
COMPUTATION 

For our computations, we use a simple model that describes the 
Galactic centre. We assume M, = 3.6 x IQ^Mq and KL = Mq 
for the mass of the MBH and the mass of each star, respectively. 
The radius of influence of the MBH, where the mass in stars is 
equal to M,, is r?i = 2 pc, and within this distance from the MBH 
there is a cusp of stars, 

iVcusp(<a) =7Vh f , (8) 

where Nh ~ Af./Af* is the number of stars within the radius of 
influence r^, and a is the slope of t he number dens ity profil e, for 
which we adopt the va lue a = 1.4 ( lAlexandeilll999l : lGenzel et al.l 
l2003l;lAlexan"de3l2005l) . 

In order to compute the torque from a cluster of stars on a 
given test star efficiently, we make use of the "wire approximation" 
suggested by Rauch & Tremaine ( 1996). We consider a test star in 
an orbit of some given initial eccentricity et and semi-major axis 
ttt — 0.01 pc in the sy-plane, with the orbit's peri-apse on the pos- 
itive s-axis. It is surrounded by a cluster of field stars whose eccen- 
tricities are drawn randomly from an isotropic distribution Niso{e) 
(see Eq. (4)) and semi-major axes from the distribution given in 
equation ([Sj. 

We truncate the semi-major axes of the field stars at 5at = 
0.05 pc, giving — 10000 stars. The orbits for the field stars start 
in a configuration similar to the test star and undergo a number of 
rotations: They are first rotated around the z-axis by an angle <j), 
drawn from a uniform distribution in [0, 2n); then rotated around 
their latus rectum by an angle 6, cosine of which is drawn from a 
uniform distribution in [— 1 , 1] ; and finally rotated around their ma- 
jor axis by an angle 7, drawn from a uniform distribution in [0, 27r). 
For a point starting from (a;, y, 0) the final coordinates (x' ,y' , z') 
are given by 

x — — y cos sin Q sin 7 

-I- y sin (^(2 cos^ </!> cos Q — cos 20) cos 7 (9) 
-I- X cos (^(cos 2(j) cos 6 + 2 sin^ (j)) ; 

y' = — y sin (j) sin 6 sin 7 

+ y cos 0(2 sin^ 0cos 9 + cos 2(f)) cos 7 (10) 
-|- X sin 0(cos 20 cos 9 — 2 cos^ 0) ; 

z' = — wfcos 6 sin 7-1-2 cos S sin sin 9 cos 7) 

(11) 

— X cos 20 sin 9 

To calculate the torques, we represent the orbits by discrete 
points that are equidistantly spaced in the mean anomaly of the 
orbit. We start by using 64 points on each orbit to calculate the 
torque. We estimate the error in our calculation by recomputing the 
torque with the points which are in the middle (in mean anomaly) 
of the points just used. If the relative difference between the two 
torques calculated, 5t = |( | tt. | — | T2 1) / (| ri | + | 1) | , is larger than 
0.01, we double the number of points and repeat the computation. 
Once the desirable tolerance is reached, we use r = (ri + T2)/2 
for torque. To limit the time spent for computation, we use at most 
65536 points per orbit. When we quadrupled this value during the 
test runs, we obtained virtually identical results. Typically, a few 
thousand points per field star were required. 
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Figure 1. The z component of the torque computed from stars with semi- 
major axes smaller than a given value, amax- The flattening of the curves 
implies that contribution from stars with semi-major axes much larger than 
the test stale's is negligible. Results for e = 0.999 are very close to e = 0.99 
case and are not shown here. 

4 RESULTS 
4.1 The torque 

To determine the eccentricity dependence of torque on an orbit, we 
made simulations with et = 0, 0.2, 0.4, 0.6, 0.8, 0.99 and 0.999. For 
each value of et, we carried out 80 simulations with TV = 10000 
stars each and averaged over the results. 

As a first result, we confirmed that the total torque is pro- 
portional to N^^'^. However, stars very far from the test star will 
not exert any discernible torque. We therefore first determine be- 
yond which point the contribution from stars become negligible. In 
Fig. [T] we plot the z component of the torque on a star of semi- 
major axis at as a function of amax, where Umax is the cut-off for 
the semi-major axes of the field stars. This figure shows that stars 
with semi-major axis larger than the test star's apo-centre distance 
Tapo = at(l -I- et), contribute very little to the net torque on the 
test star. Motivated by this, we normalize the torque by 

T = ViV(<2at)^^ . (12) 
It 

Fig.[T]already shows a strong dependence of the z component 
of the torque on eccentricity. In Fig.|2]we show this more explicitly 
by plotting the z component of the torque as a function of eccen- 
tricity. As expected, the torque vanishes for e ^ 0, and has finite 
values for e > 0. We find that the result is consistent with a linear 
growth of the torque as a function of e. The best linear fit gives 

r, = /3,(e)f = 0.25ef . (13) 

In Fig.[2l a cusp with et = 1.4 was assumed. We have performed 
another set of calculations in which a = 2, which also showed a 
linear eccentricity dependence of Tz . It can thus be concluded that 
the result in equation il3t is not strongly depended on the particular 
choice of a. 

The component of the torque perpendicular to the angular mo- 
mentum, Tx = \/ T"! 4- changes its direction but not magnitude. 
Th e resulting relaxation proc ess is called vector resonant relaxation 
bv iRauch & Tremaind ( Il996 ). We plot this component as a func- 
tion of eccentricity in Fig. [5] Here, the data are consistent with the 
torque being a quadratic function of eccentricity: 

r_i = /34e) f = 0.28 (e' + 1/2) f . (14) 
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well suited to determine this dependence, and we do not consider 
tiiis possibility here. By rotating a test wire in a cluster, we find that 
typical variations of the torque occur over angles of ~ 7r/2. Moti- 
vated by thi s, we use uj — tt/2 , smal ler than the value of lo — ir 
assumed bv lRauch & Tremaind ( 1 19961) . We note that a larger value 
for UJ leads to more effective resonant relaxation. We now discuss 
the three processes leading to reorientation of the test star's orbit 
with respect to the cluster. 

The general relativistic precession time is given by 



where 



J 



JhSO 



P = 



12r; 



-{l-e')P{a), 



(15) 



Figure 2. The z component of torque (parallel to the star's angulai' momen- 
tum) as a function of eccentricity. The large error bars give an estimation 
of the root mean square variations of the torques. These variations can be 
of order unity, so that the torque of a star with given eccentricity can vary 
considerably depending on the ambient stellar cluster; the torque for a given 
configuration is expected to lie within the large en'or bars. The small error 
bars estimate the uncertainty of the average torque of a given eccentricity, 
based on 80 different configurations of the host cluster. The average torque 
is well determined: if it were to be computed again from 80 different cluster 
configurations, it is expected to lie within the small error bars. 
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Figure 3. The components of the torque peipendicular to angular momen- 
tum, as a function of eccentricity. For clarity, we do not plot the error bars 
in this figure, but they are comparable to the ones shown in Fig.|2] 



Note that the x component of the torque vanishes for large e, and x 
and y components become equal to each other for small e. 



4.2 The resonant relaxation time 

For the translation of the torques into a relaxation time, the time- 
scale for the stars to change their orientation with respect to the 
host cluster in Eq. Q needs to be determined. Three relevant 
processes are (1) precession of the test star due to general relativity; 
(2) precession of the test star due to the extended distribution of the 
host cluster and (3) precession of the host cluster itself due to its 
own extended distribution. 

Torques are assumed to be coherent when the orbit has pre- 
cessed less than an angle ui, and to make a random walk for angles 
> u. The precise value of the coherence angle cj, which determines 
the steps of the random walk, is not clear In particular, the coher- 
ence angle may itself depend on eccentricity, ui = u}{e). Our wire 
method, which does not include the evolution of the orbits, is not 



4GM. 



(16) 



is the angular momentum of the last stable orbits for eccentric or- 
bits, and rs ~ 2GM,/c^ is the Schwarchild radius of the black 
hole. 

Since the potential is not exclusively dominated by the MBH, 
but there is a contribution of the stellar cluster as well, the orbit 
of the test star precesses. T he precession rate o f a star in a cusp 
near a MBH was derived bv llvanovet"aD l l2005l) . Here we briefly 
summarize the result for a = 3/2. Let Su}{a,e) be the change 
in the angle of the peri-centre during one orbit. The timescale for 
precession due to extended mass distribution is then given by 



tM{a,e) = —- 

2duj(a, e) 

llvanov et al showed that 



5ijj(a, e) = 4 — 
where 



3/2 



2e de 



F{e), 



(17) 



(18) 



4E 



2e 



2e 



(19) 



and E and K are complete elliptic integrals. For this result, an a = 
3/2 power-law was assumed , which simp l ifies t he equations. For 
more general expressions see llvanov et aLl( l2005l) . 

Since general relativistic precession takes place in the oppo- 
site direction, the rate at which the star's orbit precesses due to the 
combined effects of general relativity and the extended potential is 
l lHopman & Alexandeni2006ai) 



C{a,e) = 



tGK{a,e) tM(a,e) 



(20) 



For some combinations of (a, e), the precession time t2j{a, e) 
can become very large, implying that the star does not precess with 
respect to inertial space. However, for the efficiency of resonant 
relaxation, it is the orientation of the orbit with respect to the other 
stars that matters. If most of the other stars do precess, the torque on 
the test star will still fluctuate. We therefore define the precession 
time of the stellar cluster as 



tt{a.) 



Cia,e = 0.7). 



(21) 
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Figure 4. Resonant relaxation time Trr as a function of eccentricity e. for 
several examples of semi-major axis, for the parameters of the Galactic cen- 
tre. For small e, T^a — oo, because /3s (e) — > 0. The resonant relaxation 
time decreases with e, mainly because the torques t\ oc e. For large e and 
small a, Tpj^p increases again with e, because the precession time becomes 
very short due to general relativistic effects. This effect is not of importance 
for very large semi-major axes. 
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Figure 5. The resonant relaxation time Trr as a function of semi-major 
axis a, for several choices of eccentricity. Far away from the MBH, the 
Trr increases with distance, but it reaches a minimum near 0.01 pc, where 
general relativistic precession starts to dominate the precession rate. For 
high eccentricities this happens farther away from the MBH then for small 
eccentricities. 



The eccentricity e — 0.7 is the median eccentricity for an isother- 
mal eccentricity DF; approximately half of the star precess more 
rapidly than t2j{a, e = 0.7), and half of the star experience slower 
precession. 

The limiting time scale for resonant relaxation is then 



t^{a,e) = min |^f^(a, e), tj^'(a) 



(22) 



Using equation J22b in equation ([7} gives the resonant relax- 
ation time as a function of a and e. 



4.3 Resonant relaxation for a simple model of a galactic 
nucleus 

We apply our results to a simple model which may describe a galac- 
tic nucleus similar to our Galactic centre. For masses, we assume 
that M. = 3.6 X 10^ Mq and M^, = Mq. The radius of influence 
of the MBH, where the mass in stars is equal to M,, is Vh = 2 pc, 
and there is a cusp of stars, with a — 3/2 (see Eq.[8}. 

In Fig. |4]we show TB.R{a,e) as a function of e for several 
choices of a. For small eccentricities, Trr becomes very large, and 
non-resonant relaxation is much more effective in changing the an- 
gular momenta than resonant relaxation. The resonant relaxation 
time then decreases with e, but for small semi-major axes it in- 
creases again near e = 1, because general relativity causes rapid 
precession of the orbit (Eq.|15l>. 

In Fig. |5]we show rRR(a, e) as a function of a for several 
choices of e. For large a, precession is dominated by mass preces- 
sion, and Trr oc a. Closer to the MBH, at a distance of ^ 0.01 pc, 
general relativistic precession starts to dominate. This happens at 
larger a when e is large. Near the minimum, general relativistic 
precession and precession due to the extended cluster of stars can- 
cel (Eq. |20t , and resonant relaxation is limited by the precession 
rate of the ambient cluster. For yet smaller a, resonant relaxation 
becomes limited by general relativistic precession. At distances of 
~ 0.01 pc, the resonant relaxation time becomes for high eccen- 
tricities as small as a few x 10^ yr. 



5 SUMMARY AND DISCUSSION 

In this paper, we have shown that the net torque of a cluster of stars 
on a test star of eccentricity e, is proportional to e (Ea.ll3t. From 
this dependence, and the dependence of the precession time on ec- 
centricity, we determine for the first time the resonant relaxation 
time as a function of e and a (equations 171 and[22t. 

Resonant relaxation may play an important r ole in sev- 
eral phenomena near MBHs. iRauch & Tremairi3 ( Il996h and 
lRauch&Ingallsl l ll998h estimated that resonant relaxation may in- 
crease the rate of tidal disruptions of stars by the MBH by a factor 
^ 2 due to the increased rate at which stars are driven towards 
the loss-cone. It has also been suggested that resonant relaxation 
has modified the distribution of the young star cluster known as the 
"S-stars" in the Galactic c entre ( lLevinl20o3 : lHopman & Alexander! 
l2006al : |Perets et al.ll2007h . 

Resonant relaxation also plays a role in the formation of grav- 
itational wave sources. Compact remnants that spiral into MBHs 
due to the emission of gravitational waves are an important po- 
tential source of gravitational waves for the Laser Interferometer 
Space Antenna (LISA). With the exception of the Galactic cen- 
tre (Hopman, Freitag & Larson 2007), such extreme mass ratio 
gravitational wave sources are not observable until they orbit on 
very tight orbits with periods less than an hour. Since such stars 
originate from orbits relativ ely close (~ 0.01 pc) to the MBH 
jHopman & Alexandedl2005h . resonant relaxation plays an impor- 
tant role in the eve nt rate, and can lead to an incr ease of nearly an 
order of magni t ude i Hopman & Alexandej2006ah . With the excep- 
tion of lFreitad ( 120011120031) ? estimates of the event rate have relied 
on semi-analytical models which were not full y two dimen sional in 
{E, J)-space. In particular. [Hopman & Alexander! l2006al) treated 
the resonant relaxation time as averaged over eccentricities. 

The determination in this paper of rRR(a, e) as a function of 
semi-major axis and eccentricity, allows for implementation of res- 
onant relaxation in Monte C arlo codes such as those presented in 
!Freitag & B"en3 ( !200ll !2002!) . The fact that the torques depend on 
e can be of considerable importance for these results. In a com- 
panion paper (Hopman & Giirkan 2007) we use the eccentricity 
dependence derived in this paper to find the steady state angular 
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momentum distribution of stars in presence of resonant relaxation, 
and address the consequences for the processed mentioned here. 
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